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Abstract 

Effects of superconducting fluctuations on the Fulde-Ferrell (FF) state are discussed in a spatially 
isotropic three-dimensional superconductor under a magnetic field. For this system, Shimahara 
recently showed that within the phenomenological Ginzburg-Landau theory, the long-range order 
of the FF state is suppressed by the phase fluctuation of the superconducting order parameter. 
[H. Shimahara: J. Phys. Soc. Jpn. 67 (1998) 1872, Physica B 259-261 (1999) 492] In this 
letter, we investigate this instability of the FF state against superconducting fluctuations from the 
microscopic viewpoint, employing the theory developed by Nozieres and Schmitt-Rink in the BCS- 
BEC crossover field. Besides the absence of the second-order phase transition associated with the 
FF state, we show that even if the pairing interaction is weak, the shift of the chemical potential 
from the Fermi energy due to the fluctuations is crucial near the critical magnetic field of the FF 
state obtained within the mean-field theory. 

PACS numbers: 03.75. Ss, 03.75.-b, 03.70.+k 
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Introduction- Under a strong magnetic field, the Fulde-Ferrell (FFlstate is known to be 
more stable tlian tlie ordinary BCS state in singlet superconductors. 

mm 

The FF 

state has a spatially oscillating order parameter of A(r) ~ e'^^'^ '" and can exist above the 
Clogston limit where the ordinary BCS state no longer exists. 6] Since this state originates 
from the inequivalence of the Fermi surfaces of electrons with up- and down-spins under 
magnetic field, the FF state can also be induced by a molecular field in ferromagnetic 
superconductors. Indeed, the FF state in the magnetic superconductor RuSr2RECu208 
(RE=Y, Gd) was recently discussed. 

BBS 

In this letter, we discuss the effects of superconducting fluctuations on the FF state. In 
ordinary BCS superconductivity, the long-range order is possible in three-dimensional sys- 
tems, while in one- and two-dimensional cases, it is absent at finite temperatures, as proved 
by Hohenberg. jiol. Although Hohenberg's theorem does not apply to the FF state, the 
FF state is also usually believed to be stable at least in three-dimensional systems. However, 
Shimahara recently pointed out that within the phenomenological Ginzburg-Landau (GL) 
theory, even in three-dimensional systems, the long-range order of the FF state is absent 
due to the effect of the phase fluctuation of the order parameter when the system is spatially 



isotropic. 



In this letter, we present a microscopic approach to the stability of the FF 



state, focusing on the phase transition temperature Tc. We show the absence of the second- 



order FF phase transition, which is consistent with the result obtained by Shimahara. 
In addition, we point out the importance of the contribution of fluctuation to the chemical 
potential in the absence of the second-order phase transition associated with the FF state. 

The theory describing the effect of superconducting fluctuations on Tc within the Gaussian 
approximation was developed by Nozieres and Schmitt-Rink (NSR theory) 1^ to investigate 
the BCS-BEC crossover problem in strong-coupling superconductors. In this theory, Tc and 
the chemical potential /i are determined self-consistently. The latter quantity is determined 
from the equation of the number of electrons including superconducting fluctuations. When 
the pairing interaction is strong, fi deviates from the Fermi energy due to the fluctuations. 
On the other hand, this effect is weak in weak-coupling BCS superconductivity. Thus in the 
latter regime, we can actually regard fi as the Fermi energy. 

Based on the NSR theory, we can expect that if the FF state is unstable against su- 
perconducting fluctuations, the chemical potential plays an important role near the critical 
magnetic field of the FF state obtained within the mean-field theory (h^p) even in the weak- 
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coupling regime. Motivated by this expectation, we apply the NSR theory to the FF state. 
Based on this framework, we show that near h^p, the shift of the chemical potential induced 
by the superconducting fluctuations is crucial in clarifying the absence of the second-order 
FF phase transition. 

Application of the NSR theory to the FF state- The NSR theory from the viewpoint of 
the Gaussian fluctuations was discussed by Engelbrecht et al. using the functional integral 
method. 15| In this letter, we extend their theory to a three-dimensional weak-coupling s- 
wave superconductor under a magnetic field. Although this extension is straightforward, the 
resulting fluctuation effect on the FF state is found to be quite different from the ordinary 
BCS case. 

The partition function in the path integral representation is given hj Z = J V'^*'D'i/(je~^ , 
where S is the action having the form 

S = l^ drj drE ^li— + 1^ - - M^<x - f/vl>|vl>Jvl/^M/^]. (0.1) 

Here, \l/(r, r) is a Grassmann variable describing electrons. Electrons feel the external mag- 
netic field (or molecular field) h and pairing interaction U. In eq. ()0.1|) . we have neglected 
the orbital effect associated with the magnetic field. Equation (jO.lj) also assumes a spatially 
isotropic system with no lattice effect. 

We introduce the superconducting order parameter A using the Stratonovich-Hubbard 
transformation as Z = J V'^lV^^VA*VAe-^ = fVA*VAe-^^ , where S = 
/o^dr/dr[^ - ^t^-i^]. In S, = (^|(r, r), ^^(r, r)), while G is given by 



A(r,r) 
yA*(r,T) 

= Go ' + A. (0.2) 



Here, cTj (i = 1, 2, 3) are the Pauli matrices. 

The mean field solution for the order parameter Amf is obtained from 5S/\/5A* = 0, 
which gives the ordinary gap equation as Amf(i") = f/(^j(r)\l/|(r)). When we consider 
the FF state described by AmfIi") = ^mfg^''^'"''", this equation can be written as Amf = 
^J2p,wm G^{p,iuJm), where {p,iujm) is the (12)-component of the one-particle matrix 
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thermal Green function including the mean field Amf: 



^ i^m — h — £p+qpp /2 + ^MF 



(0.3) 



^ ^MF i'^m — h + e__p+qpp/2 — /i 

"ermion Matsubara 

p n n r 

frequency. At Tc with h = /ip/, the gap equation is reduced to 



Here, e±p+qpp/2 = (p ± q/2)^/2m is the kinetic energy while is the 

)Q, 0,0, 1^1 

^ _ jjy^^ ^ ^ /(gp+qFF/2 — + ^FF ) ~ /(g-p+qFF/2 ~ M ~ ^FF ) 
p £^p+qFF/2 + £^-p+qFF/2 ~ 

where f{e) is the Fermi distribution function. Within the mean field approximation, eq. 
fl0.4|) gives the well-known h — T phase diagram shown in the inset in Fig. 1. 

m 

Now we consider the fiuctuation effect around the mean field solution in the normal state 
(Amf = 0). Executing the functional integral in terms of and \E'* in the partition function 
Z, we have 

ri^ r lAP 
5a = / dr / dr^ -TrloghG-^] 



J U 

= / dr/dr^-Trlog[-Go-l]-5:-Tr[-GoA]^ (0.5) 
Jo J U ^^n 

In the following, we retain the summation in the last term up to ri = 2 (Gaussian approxi- 
mation). The first-order term (n = 1) vanishes because Amf = above Tc. Now, we write 
A(r,r) = A^^Vi — A'^^V2 in eq. ()0.5|) . where A*-^) and A^^^ are, respectively, the real and 
imaginary parts of A(r, r). When we introduce the (p, tu^) -represent at ion instead of (r, r)- 
representation, eq. ()0.5|) under the Gaussian approximation is reduced to (g = (q, ii/^), 
where Vm is the boson Matsubara frequency) 

f^u , y -f/n2i(g) l + f/n22(g) j l^Af ] 

Here, A^-'-' is the Fourier-transformed A*^-'''(r, r) and Hjj is the generalized polarization func- 
tion defined by[3l '^iM) = ^T,p TT[Go{p)(JiGo{p + q)(yj], where Go{p)~^ = ium -h-{ep- 
fi)crs- In the superconducting state, when we take the order parameter proportional to the 
cTi-component, Hn and n22 represent, respectively, the amplitude and phase fiuctuations 
of the order parameter, while ni2 and n2i describe their coupling. In the normal state, 
they satisfy Hn = n22, ni2 = -n2i, Iljj{-q) = Ujjiq) and ni2(-g) = -ni2(g). After 
the functional integration in terms of A and A* in the partition function Z, we obtain the 



thermodynamic potential Q = —TlogZ = Qq + TJ2q^og{l — UIl{q)) , where Qq is the ther- 
modynamic potential for free electrons; n(g) = — (IIii + 11122) describes superconducting 
fluctuations in the normal state: jl^ 

n^^) _ ^ 1 - /(gp+q/2 -fl + h)- /(£p+q/2 ' ^^ ' h) 
p ^P+q/2 + £^-p+q/2 — 2fi — iUm 

When we take h = 0, the thermodynamic potential Q is reduced to that obtained by Nozieres 
and Schmitt-Rink. [ij] 

The chemical potential is determined by the equation of the number of electrons, which 
is obtained from the identity = —dQ/dfi: 

"^""^"^/^V i-^n(g) a/. 

= + - / dzN-B{z)lm[- — — — : — — 

= Nl + N^i^. (0.8) 

Here, = Sp.o- fi^p — o'h) is the number of free electrons and Nb{z) is the Bose distribution 
function. In eq. (jO.Sj) . A^fl is the fluctuation contribution to N. 
For later convenience, we rewrite eq. ()0.4|) as 

l = [/n(qFF,ii^m = 0). (0.9) 

The transition temperature Tc and fi are determined self-consistently from eqs. ()().8j) and 

Stability of the FF state- In weak-coupling three-dimensional BCS superconductors, the 
fluctuation effect described by A^fl is weak; thus we can neglect A^fl- Then, since eq. ()().8|1 
is reduced to the equation of the number of free electrons, we can safely take ^ = ep, where 
Ep is the Fermi energy of free electrons ai h = 0, as long as the temperature and magnetic 
field dependence of /i can be neglected. 

However, the situation is quite different in the FF case. At h^p, we find from eq. ()0.9p 
that H(q, 0) must be maximum and be equal to 1/U at qFF- As an example, we show in Fig. 
1 the momentum dependence of H(q, 0) at h^J{T/T^ = 0.1), where T° is the phase transition 
temperature at h = 0. We find that eFn(q, 0) = 2.5, i.e., ?7H(q, 0) = 1, at |q| = 0.024A;f, 
where kp is the Fermi momentum. (For comparison, the case of /i^/(T/T° = 0.7), at 
which the BCS state appears, is also shown. In this case, H(q, 0) is maximum at q = 0.) 



When we expand n(q, ium ^ z + iS) around qpF and z = at h^J, we obtain (note that 
l-[/n(qFF,0) =0) 

1 - f/n(q, z) - C(qFF)(|q| - |qFF|)' - iaz. (0.10) 

Here, C(qFF) = -U/2 ■ d^U(QvF,0)/dq'^ > and a = {'KU/8T)N{ij)sech.^h/2T, where N{fi) 
is the density of states at /i4l7j We note that since the system is isotropic, the expansion in 
eq. ()0.10p includes (|q| — |qFF|)^, but not (q — qFF)^- Substituting eq. ()0.10p into eq. (jO.Sp . 



we obtain 



77 „ J-oo C(qFF)(|q| - IqFFi) - OfJ. 



ivg+^^"'r-°>]i:r^.- ' 



77 djJ, q J-oo C(qFF)^(|q| - IqFFi)'' + 

TU gn(qFF,0) 1 

^^c(qFF) V(lql-lqFFl)2' 



(0.11) 



In the second hne, since the region satisfying z ~ and |q| ~ |qFF| is important, we ap- 
proximately set Nb{z) T/z and (9n(q, z)/dfi ^ (9n(qFF, 0)/dfi. The last expression in eq. 
(10.1111 always diverges at finite temperatures as long as |qFF| 7^ 0, even in three-dimensional 
systems. Namely, no solution exists that satisfies both eqs. fl().8|l and Since the diver- 

gence comes from A^fl, this result indicates that the superconducting fluctuations suppress 
the second-order phase transition associated with the FF state. 

Besides the FF state, various nonuniform superconducting states which are described by 
A(r) = AJ2j^^'^^ '^ ~ const.) have been proposed under magnetic field (for example, 
A(r) = Acos(q ■ r)).]^, Q] Although these states have different free energies below the 
critical magnetic field (within the mean- field theory), their second-order phase transitions 
are still determined by eqs. ()().8|1 and ()0.9j) in spatially isotropic systems. Thus, as in the case 
of the FF state, their second-order phase transitions are also absent due to the fluctuation 
effect. 

We also note that the factor l/(|q| — |qFF|)^ in eq. (jO.llI) leading to the divergence 
originates from the fact that all the states satisfying |q| = |qFF| are degenerate in the 
present isotropic system. When we take into account a crystal lattice, the band structure 
and/or the pairing interaction may be anisotropic, which leads to the anisotropy of n(q, iz/m) 
in momentum space. Then the direction of the most stable FF state may become discrete, 
reflecting the crystal symmetry. In this case, the factor l/(|q| — |qFF|)^ in eq. ()0.1ip is 



replaced by l/(q — qpr)^ around qpF- Then, by the transformation q — qpF — ^ q, the 
stabihty of the FF state becomes similar to that in the case of the ordinary BCS state; the 
divergence is absent in three-dimensional systems. Namely, the FF state becomes stable in 
spatially anisotropic three-dimensional systems where the direction of the most stable FF 
state is discrete. 

It is worth noting that the difference between the ordinary BCS state and the FF state 
against superconducting fluctuations originates from the spatially oscillating character of 
the latter's order parameter, A(r) ~ Ae^'^^''' '". This wavy structure is similar to a "crystal of 
layers" stacking in the qpF-direction, such as a smectic liquid crystal. [3| This type of crystal 
is known to be unstable against structural fluctuations in three-dimensional systems, 
Thus a similar situation is also expected in the FF state. Indeed, assuming qpF = (0, 0, ^ff) 
and considering the fluctuations around qpF, we find (q_L = {qx, Qy), Sqz = Qz — Iff) 

TU an(qFF,o)^ 1 
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where we have dropped unimportant factors in the last expression. Equation ()().12|1 is 
identical to the expression for the effect of structural fluctuations of a layered crystal stacking 
in the z-direction in a three-dimensional system. I20I . |^ This similarity clearly shows that 
"structural fluctuations" of the wavy structure of the FF order parameter are the key to 
the absence of the long-range order of this state. We also note that this kind of fluctuation 



2l| 



effect also appears in a spiral magnetic state in magnetic superconductors. 

Fluctuation contribution to fi- Next we discuss what actually occurs near ^ using eqs. 
()().8|1 and Near h^p, when n(q, 0) is maximum at qmax (— qFF), eq. ()().8j) can be 

evaluated as 

" ^ 5/i V(l-^n(q„ax,0)) + C(q,„ax)(|q|-|qmax|)2' ^' ^ 

As the magnetic field is lowered to approach h^p, the factor 1 — [/n(qjnax, 0) becomes small, 
so that the fluctuation contribution to described by the second term in eq. (jU.13j) (A^fl) 
takes a large positive or negative value depending on the sign of 9n(qjnax, 0)/(9/i. When 
9n(qmax, 0)/9;U > 0, in order to satisfy eq. ()0.13|) . /i becomes small so as to decrease A^p. 
The decrease of also decreases n(qmax, 0) because of 9n(qmax, 0)/dfJ. > 0; thus the factor 



1 — t/n(qmax — * qFF,0) is still finite (or eq. ()U.9|) is not satisfied) at h^J and the phase 
transition to the FF state is avoided. Thus, the shift of the chemical potential due to the 
superconducting fiuctuations is crucial near h^p, in contrast to the case of the ordinary 
weak-coupling BCS state. 

The above discussion is also applicable to the case of dll{q^^ax, 0)/dfi < (negative A^fl)- 
In this case, in order to increase A^p, fi becomes large and thus n(qmax, 0) decreases. As a 
result, we obtain 1 — f/n(qmax Qff^O) > at h^p and the FF phase transition is again 
avoided. 

We note that this mechanism also works in the region h < h^p until n(q, 0) takes a 
maximum value at q = 0. Thus, the second-order phase transition into the FF state does 
not exist in the h — T phase diagram. 

Summary and discussion- In conclusion, we have discussed the superconducting fiuctua- 
tions in a three-dimensional superconductor under an applied magnetic field using the NSR 
theory. Under a strong magnetic field, the second-order phase transition associated with the 
FF state has been examined within the mean-field theory. However, this phase transition 
actually does not occur at finite temperatures due to superconducting fiuctuations unless 
the effect of the crystal lattice is taken into account. We also showed that, even if the pair- 
ing interaction is weak, the shift of the chemical potential induced by the superconducting 
fiuctuations is crucial for understanding what actually occurs near h^p . This is in contrast 
to the ordinary BCS case, where we can set ii = ep a.s far as the weak-coupling regime is 
considered. 

When the second-order phase transition associated with the FF state is absent, it is 
interesting to consider what sort of superconducting phase transition would be realized 
under the strong magnetic field. In this regard, it is known that within the mean field 
theory, when the FF state is neglected, the first-order superconducting phase transition 
occurs at the Clogston limit. In addition, the present discussion is within the second-order 
phase transition to the FF state and thus the possibility of the first-order phase transition is 
not denied. Since the present approach is valid only for T > Tc, we must extend the present 
method to the superconducting state in order to discuss these first-order phase transitions. 
This problem remains for future discussion. 

Although real superconductors are more or less anisotropic due to the presence of a crystal 
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lattice, we can expect spatially isotropic superfluidity in trapped gases of Fermi atoms such 



as and ^Li.j2a, 



23l 



2J| In these systems, two hyperflne states corresponding to electrons 



with up- and down-spins feel different potentials in a magnetic trap potential. This situation 
is similar to superconductivity under magnetic fleld. In addition, the spatial isotropy can 
be controlled experimentally by introducing an optical lattice to trapped gases, j^l Thus, if 
superfluidity is realized in this system, it is expected that we can examine how the FF state 
would be stabilized when the anisotropy of the system is introduced. 
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Fig.l Y.Ohashi 

FIG. 1: n(q, 0) as a function of |q| at the superconducting phase transition obtained within the 
mean- field theory (i.e., neglecting superconducting fluctuations). The solid line is at h^piT/T^ = 
0.1) ('A' in the inset) while the dashed line is at hfJ{T/T^ = 0.7) ('B' in the inset). Inset: h-T 
phase diagram within the mean-field theory. The upper line is the FF phase transition while the 
lower one is the BCS phase transition under the assumption of second-order phase transition. We 
set ^ = ep and U/ep = 0.4. (Thus eq. ()0.9() is satisfied at eFn(q, 0) = 2.5.) We have introduced the 
Gaussian cutoff 6"^*^^/^^^ with D = 2ep in solving eq. HU.4|) . The first-order transition (Clogston 
limit) and the phase boundary between the FF state and the BCS state are not shown. 
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